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It was pointed out by Tewari and Sau that chiral symmetry (H h-> —H ife f> h) of the Hamiltonian 
of electron-hole (e-h) excitations in an TV-mode superconducting wire is associated with a topological 
quantum number Q £ Z (symmetry class BDI). Here we show that Q = Trrhe equals the trace of 
the matrix of Andreev reflection amplitudes, providing a link with the electrical conductance G. 
We derive G = (2e 2 /h)\Q\ for \Q\ = N,N - 1, and more generally provide a Q-dependent upper 
and lower bound on G. We calculate the probability distribution P(G) for chaotic scattering, in 
the circular ensemble of random-matrix theory, to obtain the Q-dependence of weak localization 
and mesoscopic conductance fluctuations. We investigate the effects of chiral symmetry breaking by 
spin-orbit coupling of the transverse momentum (causing a class BDI-to-D crossover), in a model of 
a disordered semiconductor nanowire with induced superconductivity. For wire widths less than the 
spin-orbit coupling length, the conductance as a function of chemical potential can show a sequence 
of 2e 2 /h steps — insensitive to disorder. 

PACS numbers: 74.45. +c, 03.65.Vf, 73.23.-b, 74.25.fc 



I. INTRODUCTION 

The classification of topological states of matter, the 
socalled "ten-fold way" P has five topologically nontriv- 
ial symmetry classes in each dimensionality!^ For a one- 
dimensional wire geometry, two of these five describe a 
topological superconductor and the other three a topo- 
logical insulator. Each symmetry class has a topological 
quantum number Q that counts the number of protected 
bound states at the end of the wire. These end states 
are of particular interest in the topological superconduc- 
tors, because they are pinned at zero excitation energy by 
electron-hole symmetry and are a condensed matter re- 
alization of Majorana fermionsP Signatures of Majorana 
zero-modes have been reported in conductance measure- 
ments on InSb and In As nanowires, deposited on a su- 
perconducting substrateP^ 

A key distinction between superconducting and insu- 
lating wires is that Q £ Z2 is a parity index in a topo- 
logical superconductor, while all integer values Q £ Z 
can appear in a topological insulator. In other words, 
while there can be any number of protected end states 
in an insulating wire, pairs of Majorana zero-modes have 
no topological protection. The symmetry that prevents 
the pairwise annihilation of end states in an insulating 
wire is a socalled chiral symmetry of the Hamiltonian: 
H i->- —H upon exchange a «-» f3 of an internal degree of 
freedom, typically a sublattice index. 

In an interesting recent development P Tewari and Sau 
have argued (motivated by Ref. [5]) that an approximate 
chiral symmetry may stabilize pairs of Majorana zero- 
modes in a sufficiently narrow nanowire. The symmetry 
H h- > —H when e ■<-> h involves the exchange of electron 
and hole indices e,hP It is distinct from electron-hole 
symmetry, which involves a complex conjugation H t— > 
—H* and is a fundamental symmetry of the problem. 
The combination of chiral symmetry and electron-hole 
symmetry promotes the superconductor from symmetry 
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FIG. 1: Superconducting wire (S) connected at both ends to 
a normal metal reservoir (N). The current / flowing from the 
normal metal (at voltage V) into the grounded superconduc- 
tor gives the conductance G = I/V of the NS junction. The 
wire is assumed to be sufficiently long that there is negligible 
transmission from one end to the other. Chiral symmetry then 
produces a topologically protected quantum number Q £ Z. 
Both G = (2e 2 / h)Tr rher^ e and Q = Trrhc are determined 
by the Andreev reflection matrix r^ c of the junction. While 
the NS junctions at the two ends of the wire can have in- 
dependently varying conductances G and G' , the topological 
quantum numbers are related by Q' — — Q. 

class D to symmetry class BDI, extending the range of 
allowed values of Q from Z2 to Z. 

In this paper we investigate the consequences of chiral 
symmetry for the electrical conductance of the supercon- 
ducting nanowire, attached at the end to a normal metal 
contact. (See Fig. [I]) The conductance G is determined 
by the matrix The of Andreev reflection amplitudes (from 
e to h, at the Fermi level), 

G=^Trr ho rL, (1) 

at low bias voltages and low temperatures and assuming 
that there is no transmission from one end of the wire to 
the other end. We will show that the topological quan- 
tum number Q £ Z in the presence of chiral symmetry is 
directly related to the Andreev reflection matrix, 

Q = Trr hc . (2) 

The intimate relation between transport and topology ex- 
pressed by these two equations allows us to make specific 
predictions for the Q-dependence of G. 
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The outline of the paper is as follows. In the next 
section we derive Eq. ^ from the general scattering 



formulation of one-dimensional topological invariants 
and obtain model-independent results for the relation be- 
tween G and Q. More specific results are obtained in Sec. 
Ill using random-matrix theory,^ under the assumption 
of chaotic scattering at the NS interface. Then in Sec.|IV| 



take the form 



we numerically st udy a microscopic model of a supercon- 
ducting nanowirej 12 * 13 ! to test our analytical predictions 
in the presence of a realistic amount of chiral symmetry 
breaking. We conclude in Sec. |V| 



II. RELATION BETWEEN CONDUCTANCE 
AND TOPOLOGICAL QUANTUM NUMBER 

In a translationally invariant superconducting wire 
with chiral symmetry, the topological quantum num- 
ber Q can be extracted from the Bogoliubov-de Gennes 
Hamiltonian as a winding number in the one-dimensional 
Brillouin zone.^In order to make contact with transport 
measurements, we describe here an alternative scattering 
formulation for a finite disordered wire (adapted from 
Ref. H0|) , that expresses Q as the trace of the Andreev 
reflection matrix from one of the ends of the wire. The 
electrical conductance G can then be related to Q by an 
inequality. 



A. Trace formula for the topological quantum 
number 

The scattering problem is defined by connecting the N- 
mode superconducting wire (S) to a normal metal reser- 
voir (N). The 2N x 27V reflection matrix r(E) relates the 
incident and reflected amplitudes of electron (e) and hole 
(h) excitations at energy E. It has a block structure of 
N x N submatrices, 
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where we have also introduced a Pauli matrix t x acting 
on the electron-hole degree of freedom. 

We assume that both time-reversal symmetry and 
spin-rotation symmetry are broken, respectively, by a 
magnetic field and spin-orbit coupling. Electron-hole 
symmetry and chiral symmetry are expressed by 



t x t(-E)t x = 



"(E) (e-h symmetry), 
'(E) (chiral symmetry). 



(4) 



Taken together, the two symmetries imply that r(E) = 
r T (E) is a symmetric matrix. For spinless particles, this 
would be a time-reversal symmetry, but the true time- 
reversal symmetry also involves a spin-flip. 

In what follows we consider the reflection matrix at the 
Fermi level (E = 0). The symmetry relations Q then 



'hh 
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(5) 



These symmetries place the wire in universality class 
BDI, with topological quantum number determined^ by 
the sign of the eigenvalues of the Hermitian matrix r x r. 
This can be written as a trace if we assume that the 
wire is sufficiently long that we can neglect transmis- 
sion of quasiparticles from one end to the other. The 
reflection matrix is then unitary, rr^ = 1. The ma- 
trix product r x r is both unitary and Hermitian, with 
eigenvalues ±1. The topological quantum number Q £ 
{—N, . . . , —1, 0, 1, . . . TV} is given by the trace 



= i Tr r x r = Trr h< 



(6) 



All of this is for one end of the wire. The other end has 
its own reflection matrix r', with topological quantum 
number Q' = |Tr T x r' . Unitarity with chiral symmetry 
relates r and r' via the transmission matrix t, 



r t 
t T r> 



(r x Sf 



(r x r)(r x t) = - (r x t)(T x r'). 



(7) 



If we allow for an infinitesimally small transmission for 
all modes through the wire, so that t is invertible, this 
implies that Trr^r = — Tr T x r' , hence Q' = —Q. 

The sign of the topological quantum number at the 
two ends of the wire can be interchanged by a change of 
basis of the scattering matrix, S <-> t z St z , so the sign of 
Q by itself has no physical significance — only relative 
signs matter. 



B. Conductance inequality 

In the most general case, the Andreev reflection eigen- 
values R n € [0, 1] are defined as the eigenvalues of the 
Hermitian matrix product rb e rt e . Because of chiral sym- 
metry, the matrix rj lc is itself Hermitian, with eigenvalues 
p n G [—1,1] and R n — p\. These numbers determine the 
linear response conductance G of the NS junction, 



G 



N 

Go J]] Rn 



G = 2e 2 /h. 



(8) 



The factor of 2 in the definition of the conductance quan- 
tum Go is not due to spin (which is included in the sum 
over n), but accounts for the fact that charge is added to 
the superconductor as charge-2e Cooper pairs. 

The Andreev reflection eigenvalues R n different from 
0, 1 arc twofold degenerate (Beri degeneracy)! 14 * 15 ! The 
eigenvalues p n are not degenerate, but another pair- 
wise relation applies. Consider an eigenvalue p of rh e 
with eigenvector ip. It follows from the symmetry rela- 
tions ([5]), together with unitarity of r, that rhefhhV'* = 
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(r c hr C cip)* = -[r ee rheip)* = -pr-kh^*- So —p is also an 
eigenvalue of r\ lc , unless r^ip* — 0. This is not possible, 
again because of unitarity, if \p\ < 1. If also p ^ 0, the 
pair p, — p is distinct. 

So we see that the /o„'s different from 0, ±1 come in 
pairs ±p of opposite sign. They cannot contribute to the 
topological quantum number Q = ^2 n p n , only the p n 's 
equal to ±1 can contribute (because \Q\ of them can come 
unpaired). Since each \p n \ — 1 contributes an amount Go 
to the conductance, we arrive at the lower bound 



G/G > \Q\. 



(9) 



The upper bound for G/Gq is trivially N, the number 
of modes, but this can be sharpened UN— \Q\ is an odd 
integer. There must then be an unpaired p n = 0, leading 
to the upper bound 



G/G < min(N,N+(-l) N -W). 



(10) 



For N = 1 these inequalities imply G/Go = \Q\, but for 
N > 1 there is no one-to-one relationship between G and 

\Q\- 

Because the sign of Q does not enter, the same in- 
equalities constrain the conductances G and G' of the 
NS junctions at both ends of the wire (since Q' = — Q). 
Otherwise, the two conductances can vary independently. 

Both inequalities Ety and (10), derived here for sym- 
metry class BDI withjQI = 0,1,2, ...TV, apply as well to 
symmetry class D with Q — 0, 1 — essentially because 
the Beri degeneracy is operative there as well™ 



III. CONDUCTANCE DISTRIBUTION FOR 
CHAOTIC SCATTERING 

A statistical relation between conductance and topo- 
logical quantum number can be obtained if we consider 
an ensemble of disordered wires and ask for the Q- 
dependence of the probability distribution P{G). For 
chaotic scattering at the NS junction we can calculate the 
distribution from a circular ensemble of random-matrix 
theory. Such a calculation was performed in Ref . ^2 for 
a superconductor without chiral symmetry (symmetry 
class D). Here we follow that approach in the chiral or- 
thogonal ensemble of symmetry class BDI. 

The assumption of chaotic scattering requires a sep- 
aration of time scales Tdweii 3> ^mixing, meaning that a 
quasiparticle dwells long enough at the NS interface for 
all available scattering channels to be fully mixed. Con- 
ceptually, this can be realized by confining the particles 
near the NS interface in a ballistic quantum dotP^In the 
next section we consider a microscopic model of a disor- 
dered NS interface with comparable dwell time and mix- 
ing time, but as we will see, the conductance distributions 
from the circular ensemble are still quite representative. 



A. Distribution of Andreev reflection eigenvalues 

We start from the polar decomposition of the reflection 
matrix in class BDI, 
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(11) 



The matrix U is an N x N unitary matrix and the N x N 
matrices T, A are defined by 



M 



r = 

m—1 

1 

A=±e 



cos a m 

COS Ci r , 



>\Q\ 





i sin OL r . 



-I Sin <Xr, 





HQI 



(12a) 
C . (12b) 



The ± sign refers to the sign of Q. (For Q = the sign 
can be chosen arbitrarily.) The symbols 1„,0„ denote, 
respectively, annxn unit matrix or null matrix. We 
have defined £ = if the difference N — \Q\ is even and 
C = 1 if N-\Q\ is odd. The M = {N-\Q\-Q/2 angles 
a m are in the interval —n/2 < a r < ir/2. 

The Andreev reflection matrix r^c = {UKU^) T has 
eigenvalues p n — sina„ (n — 1,2,... M), p n — — sina„ 
(n = M + 1, M + 2, . . . 2M), p n = 1 (n = 2M + 1, 2M + 
2, . . . 2M + |Q|), and additionally p N = if N - \Q\ is 
odd — all of which is consistent with the general consid- 
erations of Sec. IHBI 

From the polar decomposition we obtain the invari- 
ant (Haar) measure p(r) = r^dr that defines the uni- 
form probability distribution in the circular ensemble, 
P(r)dp(r) = dp(r). Upon integration over the inde- 
pendent degrees of freedom in the unitary matrix U we 
obtain the distribution P(ai, «2, ■ • • cum) of the angular- 
variables. A change of variables then gives the distri- 
bution P(Ri,R 2 , . . . Rm) of the twofold degenerate An- 
dreev reflection eigenvalues R n = sin 2 a n . Details of this 
calculation are given in App. [X] The result is 



M 



p({Rn}) oc n R m c ' i/2 (i - R m ) M n ( R t - R tf- 

m—1 i<cj=l 

(13) 

The M twofold degenerate eigenvalues repel each other 
quadratically; furthermore, they are repelled with expo- 
nent \Q\ from the \Q\ eigenvalues pinned at unity. While 
the probability of finding a small reflection eigenvalue is 
enhanced for N — \Q\ even (£ = 0), the eigenvalue Rn = 
pinned at zero for N— \Q\ odd (£ = 1) produces a square 
root repulsion. 



B. 



Dependence of conductance distribution on the 
topological quantum number 



Integration over the probability distribution ( 13 1 of the 



Andreev reflection eigenvalues gives the distribution P{g) 
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of the dimensionless electrical conductance 

At 



g = G/G = \Q\+2j2^n 



(14) 



m— 1 



The first term \Q\ is the quantized contribution from the 
topologically protected eigenvalues, and the factor of two 
in front of the sum accounts for the Beri degeneracy of 
the M eigenvalues without topological protection. 



The conductance distribution is only nonzero in the 
interval 



\Q\<g<ndn(N,N+{-l) N -\ Q \), (15) 



see Sec. II B It is a trivial delta function, P(g) — S(g — 
\Q\), when \Q\ — N, N — 1. Explicit results for small 
values of N are 



N = l: P(g)=S(g-\Q\), 

f%-|0|) if 101=1,2, 



N = 2 : P(g) 



N = 3: Pig) 



N = 4 : P(g) 



(8g) 



-1/2 



if |Q| = 0, 



if |0| =2,3, 
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l{2g)V*0{2-g) if \Q\ = 0, 
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if |0|-3,4, 
if 101 = 2, 
if 101 = 1, 
2) if |Q| = 0. 



(16a) 
(16b) 

(16c) 
(16d) 



The step function 9(x) (equal to for x < and 1 for 
x > 0) is used to indicate the nontrivial upper and lower 
bounds of the conductance. (The trivial bounds < g < 
N are not indicated explicitly.) The distributions for 
N = 3, 4 are plotted in Fig. [2) 

The first two moments of the conductance can be cal- 
culated in closed form for any value of TV, 0, using known 
formulas for Selberg integrals.^ (Alternatively, one can 
directly integrate over the BDI circular ensemble, see 
App.lBl) We find 



(G/G ) = 
Var(G/G ) = 



N{N- 1) +Q 2 
2iV- 1 

A{N 2 - Q 2 )(N 2 -Q 2 -2N + l) 
(2N-l) 2 (2N + l)(2N-3) ' 



For N — > 00 at fixed \Q\, this reduces to 



(G/Go) 
Var (G/G ) 



N 1 0-1/4 
~2 ~ 4 + 2N 
1 _ Q 2 - 1/4 
4 2iV 2 



e>(iv~ 2 ), 

0(N- 3 ). 



(17) 
(18) 

(19) 
(20) 



The reduction of the average conductance below the 
classical value NGq/2 = Ne 2 /h is a weak localiza- 
tion effect, produced by the chiral symmetry in class 
BDI. (It is absent for the class-D circular ensemble 
The variance of the conductance in the large- N limit, 



VarG- 
metry. 



(e /h) , is twice as large as without chiral sym- 



A fundamental effect of chiral symmetry is that the Q- 
dependence of moments of the conductance is perturba- 
tive in 1/iV. In the class-D circular ensemble, in contrast, 
the p-th moment of the conductance is strictly indepen- 
dent of the topological quantum number for N > p, so 
topological signatures cannot be studied in perturbation 
theoryp3I 



IV. RESULTS FOR A MICROSCOPIC MODEL 

We study a model Hamiltonian of a disordered 
two-dimensional semic onductor nanowire with induced 
superconductivityJi^M' 

R= ( JpL + u( r ) - fj) Tz + V z a x T z 
V 2m e ff / 

OL 

+ -JT-{p X CFyT z - PyCx) + &0VyTy. (21) 

This Bogoliubov-de Gennes Hamiltonian contains the 
single-particle kinetic energy (p 2 +p 2 )/2m c s, electro- 
static disorder potential U(x,y), chemical potential /x, 
Zeeman energy Vz, Rashba spin-orbit coupling constant 
a so , and s-wave pairing potential Ao- The Pauli matrices 
Oi , Ti act on the spin and electron- hole degree of freedom, 
respectively. The two-dimensional wire has width W in 
the y-direction and extends along the x-direction (paral- 
lel to the Zeeman field). We define the spin-orbit cou- 



5 



1 

^0.5 



N =3 




i 




- 


100 


iei=i 


\ \Q\=2 




101=3 ! 






\ 






en 




\ 






s 




\ 






tri an 
. OU 












\ 






40 




\ 


























20 








i 







FIG. 3: Topological phase diagram of the Hamiltonian (211 
without disorder ((7 = 0) and without any chiral symmetry 
breaking (ap y a x = 0, symmetry class BDI). The colored re- 
gions give the value of the topological quantum number Q 
in the superconducting state (Ao = 8Ew), while the black 
lines separate regions with different number of propagating 
modes N in the normal state (Ao = 0). The topological 
phase boundaries are independent of _E so . The blue line is 
referred to in Fig. [4] 



FIG. 2: Probability distribution of the conductance for 
chaotic scattering in symmetry class BDI. The distributions 
are plotted from Eq. ( 16 1 for iV = 3, 4 modes and different 



values of the topological quantum number Q. Thick vertical 
lines indicate a (^-function distribution. 



pling length l so = h 2 (m e ffa so ) 1 and confinement energy 
E w = h 2 (2m eS W 2 )- 1 . 



A. Mechanisms for chiral symmetry breaking 



Electron-hole symmetry and chiral symmetry, 



t x Ht x 



—H* (e-h symmetry), 
—H (chiral symmetry), 



(22) 



together require that H is real. While the electron-hole 
symmetry is an exact symmetry of the Hamiltonian (21 1, 



the chiral symmetry is broken by the spin-orbit term p y a x 
associated with transverse motionP 

To quantify the stability of multiple zero-energy states, 
we follow Ref. [T? and make a unitary transformation 
H i — ^ hO HU = H' with U = exp(ia x T z y/l so ). The trans- 
formed Hamiltonian, 



H' 



2m eff 

h 



U 



a s , 
2L 



A* ]Tz + V Z <T X T Z + A (7yTy 



p x [cos(2y/l so )a y T z + sin(2y/l so )a z 



(23) 



no longer contains p y and breaks chiral symmetry 
through the final term oc p x a z . For W <C l so this term 



produces a splitting SE of pairs of zero-energy states of 
order (W/l so )Eg &p , with E gap (x a so the induced super- 
conducting gap. This simple estimate is an upper bound 
on the splitting, even smaller splittings have been found 
in Refs. IT8H2T)1 We typically find in our numerical simu- 



lations that SE < 0.05 E gap for W < l so - 

There are other methods to break chiral symmetry. 
An externally controllable method is to tilt the magnetic 
field so that it acquires a nonzero component in the y- 
direction, in the plane of the substrate but perpendicular 
to the axis of the nanowirePBThe orbital effect of a mag- 
netic field (Lorentz force) also breaks chiral symmetry, 
but this is expected to be small compared to the Zeeman 
effect on the spin. Subband-mixing by a disorder poten- 
tial or a position-dependent pairing term preserve chiral 
symmetry. This leaves spin-orbit coupling of transverse 
momentum as the most significant intrinsic mechanism 
for chiral symmetry breaking and we will focus on it in 
the simulations. We find a transition from symmetry 
class D (Q G Z 2 ) to class BDI (Q G Z) if W drops below 

^so ■ 

All these considerations apply to noninteracting quasi- 
particles. Interactions have the effect of restricting Q 
to Zg, so chiral symmetry can stabilize at most 8 zero- 
modes at each end of the wireP^S For N < 8 we expect 
the universal class BDI results (in particular the conduc- 
tance quantization) to be unaffected by interactions. 



B. Class BDI phase diagram 

For an infinite clean wire with exact chiral symmetry, 
Fig. [3] shows the phases with different topological quan- 
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FIG. 4: Conductance of a disordered NS junction, calcu- 
lated numerically from the model Hamiltonian (21 1. The 



chemical potential fj, is increased at constant Vz — 90 Ew , 
Ao = SEw (blue dashed line in the BDI phase diagram of 
Fig. 131, for three different values of the spin-orbit coupling 
length l so . Each curve is for a single disorder realization (of 
strength Uo = 180 Ew)' The conductance quantization at 
2, 3, 4 x 2e 2 /h is lost by chiral symmetry breaking as W be- 
comes larger than l so . 



turn number Q G Z as a function of Zeeman energy and 
chemical potential. (A similar phase diagram is given in 
Ref. EU) The phase boundaries are determined from the 
Hamiltonian (21 1 by setting ap y a x = 0, U = 0, and de- 
manding that the excitation gap vanishes. This happens 
at 



p x = 0, p y = p n = nnh/W, n 
Vi = Al + ( f ,-pl/2m eS ) 2 ) 



1,2, 



■ N, 



(24) 



with N the number of propagating modes in the normal 
state (A = 0). 

If one follows the sequence of Q, N values with increas- 
ing /i at constant Vz, one sees that \Q\ remains equal to 
N > 1 for a range of chemical potentials (fi — Tr 2 Ew) 2 < 
V%. For example, the sequence along the dashed blue 
line is {\Q\,N) = (1, 1), (2, 2), (3, 3), (4, 4), .... In view of 
the inequality ([£]), this implies a sequence of 2e 2 /h con- 
ductance steps. The first quantized conductance plateau 
emerges when the Zeeman energy exceeds the supercon- 
ducting gap (Vz > Ao). Additional plateaus form at 
fields, for which the Zeeman energy becomes larger than 
the subband splitting. More specifically, the n-th conduc- 
tance plateau appears for V| = A 2 , + £^7r 4 (n 2 — l) 2 /4 
(" 1.2.3....:. 



C. Conductance quantization 

To demonstrate the conductance quantization we at- 
tach a clean normal-metal lead at x = to the disordered 
superconducting wire. For x < we thus set A = 
and [7 = 0. The Andreev reflection matrix is calculated 



numerically by discretizing the Hamiltonian on a square 
lattice (lattice constant a = W/20). Disorder is realized 
by an electrostatic potential U(x, y) that varies randomly 
from site to site for x > 0, distributed uniformly in the 
interval (— Uq,Uq). 

The results in Fig. [4] clearly show the expected be- 
havior: For W — l so /2 the conductance increases in a 
sequence of quantized steps, insensitive to disorder, as 
long as \Q\ G {N, N — 1}. The quantization at \Q\ > 2 
is lost for W = 2l so because of chiral symmetry break- 
ing. The very first step G — 2e 2 /h is common to both 
symmetry classes D and BDI, so it persists. 



D. Conductance distribution 

For \Q\ < N — 2 there is no conductance quantiza- 
tion, but we can still search for the Q-dependence in the 
statistical distribution of the conductance. In Fig. [5] we 
show the distribution function for N = 4, \Q\ = 0,1,2, 
calculated by averaging the results of the numerical sim- 
ulation over disorder realizations. The parameters used 
are listed in the caption. The values of the Fermi energy 
(/xjv in the normal region and us in the superconducting 
region) were chosen in order to be far from boundaries 
where Q or N changes. 

We found that the conductance distributions depend 
sensitively on the disorder strength, demonstrating that 
the scattering at the NS interface is diffusive rather than 
chaotic. This is as expected, since chaotic scattering re- 
quires a confined geometry (for example, a quantum dot), 
to fully mix the scattering channels. Still, by adjusting 
the disorder strength Uq a quite good agreement could 
be obtained with the distribution from the class BDI cir- 
cular ensemble calculated in Sec. IIIIBI Since this is a 
single fit parameter for an entire distribution function, 
we find the agreement with the circular ensemble quite 
satisfactory. 



V. CONCLUSION 

In conclusion, we have developed a scattering the- 
ory for superconducting nanowires with chiral symmetry 
(symmetry class BDI), relating the electrical conductance 
G to the topological quantum number Q € Z. In a closed 
system \Q\ counts the number of Majorana zero- modes 
at the end of the wire, but in our open system these end 
states have broadened into a continuum with other non- 
topological states. Still, the value of \Q\ manifests itself 
in the conductance as a quantization G = \Q\ x 2e 2 /h 
over a range of chemical potentials (see Fig. 

More generally even when G is not quantized, the con- 
ductance distribution is sensitive to the value of \Q\, as 
we calculated in the circular ensemble of random-matrix 
theory (see Fig. [5]). Comparison with Ref. [HJ where the 
conductance distribution was calculated in the absence of 
chiral symmetry (symmetry class D with Q € Z2), shows 
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FIG. 5: Blue histograms: probability distribution of the 
conductance of the NS junction, calculated from the model 
Hamiltonian (21 1 in an ensemble of disorder realizations. 
Each panel has the same number of modes N = 4 in the 
normal region and a different topological quantum number 
\Q\ — 0, 1, 2 in the superconductor. The black curves are the 
corresponding distributions in the class BDI circular ensem- 
ble, given by Eq. ( |16d| |. Each panel has the same value of 
l so = 2W and Ao — 8Ew- The other energy scales (in units 
of Ew) are as follows: Q = 0: /ijv = fis = 64, Vz = 14, 
U = 180; \Q\ = 1: fi N = 64, Ms = 88, Vz = 14, U = 180; 
|Q| = 2: p, N = MS = 64, Vz = 34, U = 140. 



that chiral symmetry manifests itself even when \Q\ < 1 
— so even if there is not more than a single Majorana 
zero-mode. 

The chiral symmetry is an approximate symmetry (un- 
like the fundamental electron-hole symmetry), requiring 
in particular a wire width W below the spin-orbit cou- 
pling length Z so . Our model calculations in Fig. [4] show 
that chiral symmetry is lost for W > 2l so and well pre- 
served for W < i so /2. Existing experiments^^ n InAs 
and InSb nanowires typically have l so ~ 200 nm and 



W ~ 100 nm. These are therefore in the chiral regime 
and can support more than a single zero-mode at each 
end, once the Zeeman energy becomes comparable to the 
subband spacing. 
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Appendix A: Calculation of the Andreev reflection 
eigenvalue distribution in the BDI circular ensemble 

In this Appendix we derive the probability distribution 
P({i? m }) of the Beri degenerate Andreev reflection eigen- 
values R±, . . . , Rm m the circular ensemble of symme- 
try class BDI (circular chiral orthogonal ensemble). The 
calculation follows the standard procedure of random- 
matrix theoryP^l and is technically similar to the calcu- 
lation for symmetry class D (circular real ensemble) pre- 
sented in Ref . ITT1 

The probability distribution P({R m }) is determined 
by the invariant (Haar) measure dp,(r) — r^dr = Sr, 
which for a given topological quantum number Q char- 
acterizes the uniform distribution of scattering matrices 
in the circular ensemble subject to the symmetry con- 
straints of Eq. ([5]). Since any scattering matrix in the en- 
semble can be decomposed according to Eq. |IT| ) , i.e. pa- 
rameterized in terms of the angles a m , we can transform 
the invariant measure into dfi(r) — J dpi da m . 
The piS denote the degrees of freedom of the matrix of 
eigenvectors U and J is the Jacobian of the transforma- 
tion. From this expression the distribution of the angles 
a m follows via integration over the p^s. Up to a normal- 
ization constant we have 



P({a m }) oc J Jjjdpi 



(Al) 



The polar decomposition in Eq. ( 11 ) is not unique. As 



in Ref. [TT] the redundant degrees of freedom can be re- 
moved by restricting the independent parameters pi in 
the matrix of eigenvectors U. The total number of de- 
grees of freedom furthermore depends on N as well as on 
Q. This is best seen if one considers the reflection matrix 
f in a basis where it is a real orthogonal and symmetric 
matrix, of the form 



O 









— IjV-Q 



o 1 



(A2) 



with O a 2A^ x 2N real orthogonal matrix. In this ba- 
sis the topological quantum number is given by Q — 



| Tr f . The upper-left and lower-right blocks do not 
change under an additional orthogonal transformation 
O'n+q ® O'n-q- Group division readily gives the total 
number of degrees of freedom: dim 0(2N) — dim 0(iV + 
Q)-dimO(N-Q) = N 2 -Q 2 . Since there are M angular 
parameters a m , there must be N 2 — Q 2 — M independent 
degrees of freedom pi in the matrix of eigenvectors U. 

In order to obtain the probability distribution from 
Eq. (Al) we need the Jacobian J. It can be determined 



from the metric tensor g M „, which can be extracted from 
the trace Tr SrSr^ , when it is expressed in terms of the 
infinitesimals da m and dpi (collectively denoted as dx^): 



Tr SrSr^ = Y, g^dx^dx^ . 



(A3) 



In view of the polar decomposition (11 1 one has 

WUrW* = 5WL + dL + L5W T , (A4) 



where we abbreviated 



W 



L 



r A 
a t r 



(A5) 



Unitarity ensures = d(WlI) = dWU + WdU SW = 
-SU. Substitution of Eqs. |aI| ) and ([A5]) into Tr SrSr^ = 
Tr drdr^ = Tr (W^ drW*W*~dr*W) gives 

Tr 5r5r j = Tr (dLdrf - 2LSW T L ] 8W - 2SW 2 ) . (A6) 



From the block structure of W and L we find Tr SW 2 = 
2 Tr SU 2 and 

Tr (L5W T LUW) = 2 Tr {TSUTSU T - A5UASU), (A7) 

where we have used r T = T* = T and A* = A. 

It is convenient to express Tr SrSr^ in terms of the form 
TrAA^ = \Aij\ 2 . Using rr + AA = ijv wc find 

Tr SrSr^ = TrdLdl) + 2Tr (TSU T + SUT)(TSU T + SUT)^ 
+ 2Tr (ASU - SUA)(A6U - SUA)^ 
= Ti +T 2 +T 3 . (A8) 

The first trace simply evaluates to 

N M 



The remaining two traces T 2 and T3 need to be calcu- 
lated using the block structure of T and A in Eq. ( 12 1. We 



work out the calculation for Q — 0, N even (=> £ — 0). 
The two matrices A and T are in this case fully described 
by M — N/2 blocks of 2 x 2 matrices. The two traces 
evaluate to 



\T 2 = ^2cos 2 a k {\SU 2K2k \ 2 + \SU 2k ^ 2k ^\ 2 + 2[lm(SU 2k ^ 2k )] 2 } 
fc=i 

M 

+ Y {( cos2 «fc + cos 2 a/)(|(5C/ 2fc!2i | 2 + \5U 2k ,2i~i\ 2 + \SU 2k -i.2i\ 2 + |<SC/ 2fc _i, 2/ _i| 2 ) 

k<l=l 

~2cosaicosa k Re{5Ul t2k + 5Ul L2k _ l +5Ul_ iak + 5Ul_ iak _ l )} , (A10) 

M 

|T 3 = Y sin 2 a k [\SU 2k . 2k - SU 2k ^ 2k ^\ 2 + 4 [Im {5U 2k -i,2k)] 2 } 
k=l 

M 

+ Y {( sin2 «fe + sm 2 ai){\SU 2ka i\ 2 + \SU 2k ,2l-i \ 2 + \SU 2k ^ h2 i\ 2 + \6U 2k -i,2l-i\ 2 ) 
k <i=i 

+ 4sma k smaiRe(SU 2ki2 i- 1 5U 2k ^ 12l - SU 2 k,2i SU 2k-i,2l-i)} ■ ( Al1 ) 

Like T and A the elements of the matrix SU can be grouped into separate 2x2 blocks, denoted by the block indices 
fc, I = 1, . . . , M. We first consider the block-off-diagonal part for which we can choose as independent parameters 

SU 2k:2 i, 5U 2kj2 i-\, SU 2k -i^ 2 i, SU2 k -i,2i-i, 

with 1 < k < I < M. The real and imaginary parts, denoted by 8U n , SU 1 , produce a total of 4M(M — 1) independent 
parameters. Note that SW = —SU immediately implies SU^ 2l = —SU^i 2k i &U\ k 2l = SU^ 2k , and so on. For given 
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values of k and / the contribution to Tr Sr5r^ has the form 

o[(w£, 2 2 +(^2U-i^ 

+ 2c[5[/H )2 i-1^^2fe-l,2i + &U\ k 2 i-i5U\ k _ l 2 i — ^2fe,2Z^2fc-l,2i-l — ^^2fe,2Z^^2fe-l,2/-l] i 

where we abbreviated a = 2 ( 1 — cos cos a{) , b = 2(1 + cos cos a; ) , and c = 2 sin a k sin a; . 

I 



The contribution to the metric tensor is a block matrix 
fa 



c ON 

-c a 

a c 

V c a/ 



b -c ON 

-c & 

6c 

c &/ 



where the first and the second block correspond to the 
real and imaginary parts, respectively. The determinant 
of this block matrix is 256 (sin 2 a k — sin 2 a;) 4 . This gives 
us the contribution to the Jacobian from the off-diagonal 
matrix elements 



M 



^off-diagonal = J\ ~ ^ 



(A12) 



k<l=l 



Next we consider the diagonal 2x2 blocks. Anti- 
Hermiticity of SU implies 5U$ = (i = 1, . . . , N). We 
choose the 2>M independent parameters 



^2fe-l,2fe-l) 



SUlk-i 



2 k 



The contribution to Tr 8r8r^ for a given value k has the 
form 

v [{SUl kak ) 2 + («5C/ 2fe _i,2 fe -i) 2 ] 

where v = 1 + cos 2 and w = sin 2 afc . Note that 
8U\ k _ l 2k is fully decoupled. The contribution to the 
metric tensor is 



Integration over the N 2 — \Q\ 2 — M ancillary degrees of 
freedom of the matrix of eigenvectors U only gives rise to 
an overall constant. A transformation of variables from 
a m to R m = sin 2 a m gives the distribution |j"3| ) of the 
twofold degenerate Andreev reflection values in the case 
Q = 0, N even (C = 0). The cases Q 7^ and/or N odd 
are worked out similarly. 



Appendix B: Average conductance in the BDI 
circular ensemble 



In the circular ensemble of Sec. |III B] the 2N x 2N re- 
flection matrix r is uniformly distributed in the unitary 
group, subject to the restrictions of electron-hole sym- 
metry and chiral symmetry. The average conductance 
can be calculated directly by integration over the uni- 
tary group. We give this calculation here, as a check on 
the result (17) derived by going through the distribution 



of Andreev reflection eigenvalues. 

Unitarity (rr^ = 1) implies that the expression for 
the conductance can be written equivalently as 



G — | Go Tr (l + r hc r nc + r eh r eh - r cc rj c - r hh r hh ) 



jG Tr (l - T z rr z 



(Bl) 



v 

-w 



-w 

V 



with a determinant of 4 cos 2 a k - This leads to a contribu- 
tion to the Jacobian from the diagonal matrix elements 



M 



•/diagonal = \\ (1 Sm 2 ttfc) 



(A13) 



fc=i 



The total number of independent parameters that we 
have accounted for is AM 2 = N 2 (including the M an- 
gular parameters a m ). This is exactly the number we 
expect for N even and Q = 0. Collecting all the terms 



that contribute to the Jacobian in Eq. (Al), we obtain 
the probability distribution 



Electron- hole symmetry (r — r x r*r x ) and chiral symme- 
try (r = r T ) constrain r to the form 



J / 4 OVnO T e lT ^ /4 



(B2) 



The matrix O is real orthogonal (OO t — 1). The di- 
agonal matrix T>q has entries ±1 on the diagonal with 
TtD q = 2Q, consistent with Eq. ([6]). Substitution into 



Eq. (Bl I gives 



G= iG Tr(l 



TyOV Q T TyOV Q T ) 



(B3) 



M 



A I 



P{a k ) (x Y[(l - sin 2 a k ) JJ (sin 2 a k - sin 2 a;) 2 



k=l 



k<l=l 



(A14) 



In the circular ensemble the matrix O is uniformly dis- 
tributed with respect to the Haar measure for 2N x 2N 
orthogonal matrices. The average of a product of four 
orthogonal matrices equals^ 
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(O aa Oj3b0 lc Osd) 



2N + 1 



2N(2N - l)(2JV + 2) 
1 



2N(2N - l)(2iV + 2) 



(SapSabSjsScd + <5 SpsSbd + SasSadSjS^Sbc) 

(Sapdacdjsdbd + SafidadS-fSSbc + S al S a b6 /38& cd + ^a^^ad^ fi8&bc 



+ SasSabS/SfScd + 5 a s8 ac 6 /3-fSbd) ■ 



(B4) 



The average of Eq. (Bl I becomes 



which is just Eq. (17). 



(G) - \G (<2N + 



4Q 2 - 2N 
2N- 1 



(B5) 
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